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SUMMARY 

Theoretical performance of a new space radiator concept, the gas particle 
radiator (GPR), Is calculated. The GPR uses a gas containing emitting, sub- 
micron particles as the radiating media. A transparent window contains the 
gas particle mixture around the solid radiator emitting surface. A major 
advantage of the GPR Is that large emlsslvlty (e^ > 0.8) Is achieved without 
the use of emissive coatings. A radiation heat transfer analysis shows that 

_4 

for a modest volume fraction (~10 ) of submicron particles and gas thickness 

(~1 cm) the emlsslvlty, c , Is limited by the window transmittance. Besides 
determining the emlsslvlty, the window Is the critical element for making It 
possible for the GPR to have lower mass than a tube type radiator. The window 
acts as a "bumper" to provide meteoroid protection for the radiator wall. The 
GPR was compared to a proposed titanium wall, potassium heat pipe radiator. 

For both radiators operating at a power level of 1.01 MW at 775 K It was cal- 
culated that the GPR mass was 31 percent lower than the heat pipe radiator. 

I. INTRODUCTION 

High specific power (power radiated/radiator mass), small area and long 
lifetime are the desirable characteristics of a space radiator. These charac- 
teristics will be attained If a low mass and high emlsslvlty, , that Is 
stable for long periods (7 to 10 yr) can be achieved. 

For a tube type radiator (either a heat pipe or a pumped loop) high emls- 
slvlty (e > 0.8) Is achieved by the use of emissive coatings. Adhesion and 
emissive stability of these coatings must be obtained for long periods of time 
If a tube type radiator Is to be a successful space radiator. Generally, the 
largest mass portion of a tube radiator Is the armor that must be used to pro- 
tect against meteoroid penetration. 

The gas particle radiator (GPR) Is a new concept that has the potential 
for long lifetime, high emlsslvlty with lower mass than tube radiators. 

Figure 1 Is a conceptual drawing of the GPR. A gas which contains a suspension 
of fine particles Is contained In a sealed volume between the tube radiator 





and an outer window that are separated by a distance, D. On start-up of the 
radiator a temperature gradient will exist across the gas. This temperature 
gradient will Induce a gas flow that will distribute the particles throughout 
the gas. However, this will have to be demonstrated for a successful GPR. In 
the micro-gravity of space the particles should remain In suspension. Uniform 
particle distribution Is critical to obtaining high emlsslvlty. Therefore, It 
will be necessary to prove uniform distribution can be established before a 
GPR Is a viable space radiator. If the window Is transparent to the emitted 
radiation the gas particle mixture will yield a high, stable emlsslvlty. Past 
Investigations (ref. 1 ) of a soot-containing flame, that Is similar to the gas 
particle mixture of the GPR, have yielded large emlttances. Obtaining high 
emlsslvlty without the use of emissive coatings Is a major advantage of the 
GPR. 

It would appear that the addition of the gas-particle mixture plus a sur- 
rounding window will result In a larger mass for the GPR than a tube radiator. 
However, the addition of the window "bumper" will provide meteoroid protection 
for the emitting tube. This window "bumper" means the combined thickness 
(window *- tube) may be significantly less than the tube-only thickness and 
still provide the same meteoroid protection (refs. 2 and 3). Therefore, the 
GPR may have a lower mass than a tube radiator. 

The major problem for the GPR Is finding a suitable window material. As 
will be shown, It Is the window transmittance that limits the emlsslvlty. 

There are suitable choices for the gas and particle materials. However, there 
Is a need to experimentally verify the emlsslvltles predicted In this study. 

Suitable window materials and a verified emlsslvlty are the critical 
Issues for the GPR. However, even If these Issues can be resolved satisfac- 
torily there are other design problems that will also have to be solved. It 
will be necessary to provide a seal between the window and the radiating sur- 
face. Also, provision will have to be made for the difference In thermal 
expansion of the window material and the metallic radiating surface material. 
The analysis and results for the emlsslvlty of the GPR are presented In the 
following section. A discussion of possible window materials Is presented In 
section III, and a mass comparison Is made between the GPR and tube type radi- 
ators In section IV. Finally, conclusions are presented In section V. 

II. EMISSIVITY OF GAS PARTICLE RADIATOR 

A number of studies (summarized In ref. 1) of the emlttance of soot- 
containing gases have been carried out. These Investigations yield large 
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emlttances for modest amounts of soot concentration. Seeding with small par- 
ticles has also been proposed to Increase gas absorption of Incident radiation 
(ref. 4) or as a means of shielding a surface from Incident radiation (refs. 5 
and 6). The success of the 6PR depends on large emlttance (or absorption) of 
the gas-particle mixture, similar to these concepts. In reference 4 large 
absorption was measured as a function of wavelength, x, In the range 0.2 < X 
< 1 ym for carbon, aluminum oxide, hafnium carbide, and tungsten particles 
suspended In water. Particle diameters ranged from 0.02 to 2 ym. 

The total emlsslvlty, e , for the GPR Is derived In Appendix A. The 
most Important approximations used In this derivation are that the gas, parti- 
cles, and radiator surface are at the same temperature and the absorption 
coefficient, k^, of the gas-particle medium Is uniform and Is determined by 
the particle optical properties. It Is assumed that the gas Is transparent to 
the emitted radiation. Also, the emlsslvlty of the radiator surface, c p . Is 
neglected. Two ranges of particle size are considered In Appendix A. For 
particles of small radius, r^, (r^/X < 1) the following absorption coefficient 
Is assumed, 



(1) 


where K Is a constant, X Is wavelength and tp Is the volume fraction of 
particles In the gas-particle medium. Equation (1) Is an absorption coeffi- 
cient that fits the emlttance results for soot-containing gases (refs. 1 and 6). 
For large particles (r^/X > 1) the absorption cross section, for the parti- 
cles Is assumed to be equal to the particle cross-sectional area, irr^. There- 
fore, the absorption coefficient In this case Is the following. 

- vH - 1 f: - r o A » 1 (2> 


Where n,j Is the particle number density. 



(3) 


For the emission temperature range of Interest for the GPR (300 to 1200 K), 
more than 90 percent of the radiation will be In the range 1 < X < 75 ym. 
Therefore, the small particles must satisfy r^ < 1 ym and the large particles 
are defined by r^ > 75 ym. 
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As mentioned above, It was assumed that the gas and particles are at the 
same temperature. The validity of this approximation Is examined In Appendix B. 
For small particles (r^ < 8,^ where 8. g Is the mean free path) the heat trans- 
fer between the gas and particles will be determined by free-molecular flow, 
whereas for large particles (r^ > 8, g ) the heat transfer Is determined by con- 
tinuum flow. Results for the temperature difference, AT = T - T ,, where 

g rad 

T ra( j Is the particle temperature at a gas pressure, p^ = 0.1 atm, and a par- 
ticle temperature, T . = 1000 K are calculated In Appendix B. For both small 

rad 

particles (free molecular flow) and large particles (continuum flow) AT < 0.01 

Trad' Under these conditions of pressure and temperature It appears reasona- 
ble to neglect the difference In temperature between the gas and the particles. 

An Important part of the e T calculation Is the window specular trans- 
w w 

mission, t. . For most window materials t.. can be approximated by a square 

^ ”■ w 

wave (ref. 7). Therefore, equation (A31) was used to approximate t^. 

Using this window approximation and the assumptions described above the 

following results for Cj were derived In Appendix A for large particles, 

•t - - e "“) F o-X T • r d > 75 (4) 
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( 6 ) 
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The quantities F rt T and Q[x |w] are Integral functions given by 
0- *-u 1 p 

equations (A36) and (A40). Other quantities appearing In equations (4) to (7) 

are the speed of light, c, Plank's constant, h, and Boltzmann's constant, k, 

w 

(hc/k = 14,388 ym-K) . The window transmission, -r w » Is the magnitude of -c\ 
between the lower, x^, and upper, x , cutoff wavelengths on the window trans- 
mission, (eq. ( A31 ) ) , and Tp Is the temperature of the emitting radiation 
surface. 
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In figure 2, cj/t w for large particles (eq. (4)) Is shown as a function 
of \ u T p for several values of B. For large values of B (>4) the exponential 
term In equation (4) Is negligible and e T ■* t F , _ . In this case e T Is 

T w °-v P T 

determined entirely by the window transmission. For small particles (eq. (5)) 

the quantity c_/t Is shown as a function of XT for several values 
T W Up 

of y- Similar to the case for large particles, when y Is large (y > 1) 

e- -» t F , _ . Therefore, for both small and large particles the total 
T W 0-\uTp 

emlsslvlty will be determined by the window transmission If the parameters B 
and y are large. 

Consider the volume fraction, <p, necessary to attain large B (>4) and 
large y (>1). In the case of large particles (eq. (6)) the smallest r d 
yields the largest value of B. Therefore, assume r^ « 75 ym and for a prac- 
tical GPR, D « 1 cm. As a result <p > 0.01 Is required to attain B > 4. Now 
consider small articles (eq. (7)). The parameter, K, Is In the range 4 to 
6 ym/ym for soot (ref. 1). Therefore, assume K * 5 ym/ym. Also, considering 

the lowest temperature of Interest T * 300 K and again using D « 1 cm the 

P _a 

volume fraction necessary to attain y > 1 Is «p > 2.7x10 . A much smaller 
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volume fraction (<p > 10 ) Is required to attain y > 1 lor small particles 

than to attain B > 4 for large particles ((p > 0.01). The smallest <p Is 
most desirable for the GPR since the total particle mass will be a minimum. 
Therefore, the use of small particles (r^ < 1 ym) will yield the same with 

less mass than with large particles (r. > 75 ym) . 

d -4 

Based on the above discussion, <p > 10 , should Insure that the maximum 

feasible emlsslvlty /e T ■* t F , T \ will be attained for small particles. 

I T w °-Vp) 

However, as figures 2 and 3 show, only for X y Tp > 7X10" 1 ym K will e j/ T w be 
large (e T /x > 0.8). At T « 300 K this means X > 20 ym. There are few 
materials with good transmission at wavelengths this long. Therefore, the GPR 
will be more appropriate for higher temperatures. Window materials that can 
possibly be used for the GPR will be discussed In the next section. 

III. WINDOW MATERIALS 

As. shown In the previous section It Is the window transmission that 
determines the GPR emlsslvlty. Therefore, It Is critical to have a window 
with large transmittance In the wavelength range of Interest. Since the radi- 
ator temperatures of Interest are 300 to 1200 K the majority of emitted radia- 
tion will be In the Infrared (X > 1 ym) . Therefore, window materials such as 
ordinary glass will not be suitable. 
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Several materials are listed In Table I that have excellent transmittance 
In the Infrared. The optical transmlttances of all materials except the heavy 
metal fluoride glasses, (HMFG), were obtained from reference 7. All the other 
properties except the yield stress were obtained from reference 8. The yield 
stresses were obtained from the noted references. All properties of the HMFG 
systems were obtained from reference 9. Most of the materials shown In Table I 
are alkali halides, which are soluble In water. Although water solubility 
would be a serious problem for operation In a water-containing atmosphere, 
operation In the vacuum of space should relieve this problem. Manufacturing 
large pieces of the alkali halides of the required optical properties may be a 
problem. However, continuous single crystal fibers of cesium bromide and 
other alkali halides with excellent transmittance have been manufactured 
(ref. 10). 

The window material In the GPR must be strong enough to contain the gas 
particle mixture. Since very thin windows (<1 mm) are desirable the yield 
stress of the material will determine the limiting thickness. Yield stresses 
of the pure alkali halides Is low. However the addition of halide Impuri- 
ties (refs. 11 and 12) results In significant Increases In the yield stress. 
Whether the optical properties are altered by these Impurities Is not stated 
In references 11 and 12. 

If the GPR Is to become a workable concept an appropriate window material 
must be developed. Alkali halides are possible candidates. However, manufac- 
turing thin sheets of these materials with sufficient yield strength must be 
demonstrated. 

IV. MASS COMPARISON BETWEEN GAS PARTICLE RADIATOR AND TUBE RADIATOR 

In an earlier section the conditions necessary to obtain large emlsslvlty 
for the GPR were described. Besides large emlsslvlty, low mass Is the other 
Important characteristic of a space radiator. In this section a mass compari- 
son Is presented between the GPR and tube type radiators. 

The total mass, My of the GPR can be written as follows, 

M t = « t A t ♦ M a (8) 

2 2 
where ay (kg/m ) Is the specific mass of the radiating surface. Ay (m ) Is the 

total radiating surface area and M Is the mass of auxiliary components such as 
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the radiator working fluid and manifolding. Similarly, the total mass of a 
tube type radiator can be written as follows. 


M = a M + M 
o oo a 


(9) 


Dividing equation (8) by (9) the ratio of total masses Is obtained. 

\ 


a T A T 

M ft = <*A 
0 0 0 
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°T 
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(10) 
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°/ 


The auxiliary masses (M and H ) are assumed proportional to the radl- 

a 3 

atlng areas. 


M = A T 
a T T 


M = a 1 A 
a « 0 T 
o 


(11a) 

(lib) 


Assume the GPR Is formed by adding a gas particle mixture surrounded by a win- 
dow to the tube type radiator. In that case a relation for the ratio of spe- 
cific masses, a^./a Q , Is derived In Appendix C. The result Is given by equa- 
tion (C14). The area ratio, A^/A^, Is determined by the power radiated. For 
a space radiator redundant area must be Included to make up for radiator area 
that Is lost as a result of meteoroid penetration. As mentioned In the Intro- 
duction, the window for the GPR acts as a bumper to protect the radiator sur- 
face. Therefore, If the window Is penetrated by a meteoroid the radiating gas 
particle mixture Is lost but the radiator surface will remain operable. This 
remaining radiator surface will have an emlsslvlty, e|, that Is less than the 
the GPR emlsslvlty, e^. For the tube type radiator all area that Is affected 
by a meteoroid penetration Is lost as a radiating surface. The radiator area, 
A $ , with emlsslvlty, e Q , that remains operative at the end of the mission must 
be sufficient to emit the entire radiator load. In the case of the GPR the 
radiator load at the end of the mission Is shared by two areas. They are the 
area, A $ with emlsslvlty c^. that does not suffer meteoroid penetration and 
the area, A T - A $ , with emlsslvlty c| that has lost the gas particle mixture 
as a result of meteoroid penetration. If both the GPR and the tube type 


7 



radiator are operating with the same load and with the same view factor the 
following relation Is obtained for the area ratio, A-j./A q . 



( 12 ) 


In obtaining equation (12) the same redundancy factor, r, was assumed for the 
GPR and the tube type radiator. 


r 



Using equations (C14) and (12) the following result Is obtained. 



(13) 


(14) 


where t Is the thickness of the radiator tube of the GPR, t Is the thlck- 
p o 

ness of the radiator tube of the tube type radiator, t Is the thickness of 

w 

the window, t Is the thickness of the window attachment member, and y , y , 
s w g 

and v s are defined In Appendix C. 

The y (D/t ) term Is the mass of the gas particle mixture, the y 
g P w 

(t /t ) term Is the mass of the window, and the y (t /t ) term Is the mass 
w p s s p 

of the attachment member. The y term can be neglected In equation (14) for 

-4 ^ 

small volume fractions (<p < 10 ). Referring to equation (C17), y (D/t ) ~ <p 


(D/tJ, since P « P and 
-2 y 


p 6 * V 


Therefore, If <p « 10 


-4 


and D/t p * 10\ 


v g (D/t p > « 10 

As mentioned In the Introduction, the window of the GPR will act as a 

bumper against meteoroid penetration of the radiating surface. Therefore, the 

total thickness, t + t^, can be significantly less than t Q and still provide 

the same meteoroid protection. For high speed meteoroids (~20 km/sec) the 

bumper fragments and vaporizes the meteoroid dispersing the fragments over a 

large area so that penetration of the main wall does not occur. Data from 

Explorer 46 (ref. 3) using a stainless steel bumper and main wall Indicated 

that t + t « t /6.9 provided the same meteoroid protection as t . In this 
p w o K o 

case t w Is the bumper thickness. Since meteoroid penetration depth Is not 
sensitive to the target material (refs. 13 and 14) It Is expected that similar 
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behavior will exist for other materials. In reference 3 the optimum distribu- 
tion of mass between the bumper and the main wall was calculated to be 

t /t = 0.1 to 0.2. Using this result together with t + t = t /6.9 the 
w p p w o 

following Is obtained. 

t 

0.12 < ^ < 0.13 (15) 

o 

Therefore, unless e /e„. » 1, the magnitude of ah /a A (eq. (14)) 

0 T T T 0 0 

will be small If t /t Is given by equation (15). If a T A T /a A Is small 

po TToo 

then a significant reduction In the mass ratio M/M , (eq. (10)) will result. 

T o 

Now consider a comparison between a proposed heat pipe radiator (ref. 14) 

and a hypothetical GPR. A cross section of this radiator and this radiator 

with a gas-particle mixture and window added Is shown In figure 4. The heat 

pipe radiator of reference 14 uses a flat plate radiating surface of titanium 

at Tp = 775 K with an emlsslvlty, c Q = 0.9 and a total radiated power of 

1.01 MW. The radiator wall thickness (t = 0.60 mm) and redundancy, 

o 

(r = 321/360) was designed to provide a lifetime of 7 yr with 0.99 probability 
of survival against meteoroid penetration. 

The emitting gas particle mixture for the hypothetical GPR consists of 

3 

small carbon particles (p^ = 1.6 gm/cm ) and helium. Helium will be trans- 
parent to the emitted particle radiation. Also, Its low molecular weight and 
thus high heat transfer coefficient (eq. ( B1 9 ) ) means there will be a neg- 
ligible temperature difference between the helium and carbon particles. Since 
It Is desirable to make the window as thin as possible the helium gas pressure 
should be as low as possible. The helium pressure was set to obtain 

AT < 0.01 T ., where AT = T - T . Is the temperature difference between 
rad g rad 

the helium gas (T ) and the particles (T .). Using equations (B19), (B27a), 

9 r 3Q 

and (B28) for T = 775 K and e /a * 1 , the helium pressure requirement Is, 

2 9 s 

p g > 2000 n/nr (15 torr). 

Window material was chosen to obtain high emlsslvlty. From figure 3 for 

4 

Y > 1 It can be seen that \ u T p >10 ym K for c t /t w > 0.9. Therefore, for 

T = 775 K, \ >13 ym. Referring to Table I we see there are several mate- 
p u “ 

rials that satisfy > 13 ym and also have a melting point significantly 
above 775 K. For the hypothetical GPR sodium chloride was chosen as the window 
material. Whether NaCl can be produced In large area, thin sheets suitable for 
a GPR Is a question that Is yet to be answered. 

Consider the emlsslvlty that Is possible with the carbon particle, helium 
mixture, The parameters, T (eq. <7» , and XT p determine 
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the emlsslvlty, c t /t w . For D = 1 cm, K = 5 ym/ym and T p = 775 K, then y > 1 

for <p > 1.03xl0“ 4 . Therefore, as long as <p > 1.03xl0 -4 the y > 1 curve In 

figure 3 can be used to obtain c t /t . For NaCl X = 15 ym and t > 0.9 

I w » u w 

(Table I). Using figure 3 for X u T p = 1.16x10^ ym K, y > 1, and t w = 0.9 

results In = 0.84. Thus c^. Is slightly less than the assumed emlsslvlty, 

c q = 0.90, of the heat pipe radiator. When the area ratio, A^./A q (eq. (12)), 

Is calculated using c| = 0.3 and r = 321/360 however, It Is found that 

A t /A = 1.03. The Increase In area resulting from e /e_ >1 Is nearly 
To o T Jj 

compensated for by the redundant area savings, [1 + ej/e-j. (1/r - 1)] term 
In equation (12) . 

As already mentioned the thickness of the heat pipe wall Is t Q = 0.6 mm. 

Applying equation (15, results In 0.072 < t p < 0.078 .. Although It may be 

possible to use such a thin wall thickness, t = 0.1 mm was arbitrarily chosen 

as a practical limit on t . This same limit was also chosen for t , so that 

p w 

t = t = 1/6 t = 0.1 mm. With t + t = 0.2 mm meteoroid protection will 
w p o w p 

be more than required, (t + t » t /6.9). Titanium of thickness, t = t = 

M w p o s p 

0.1 mm, will also be assumed for the connecting member between the plate and 
window. 

As well as providing meteoroid protection, the window and emitting plate 

must be of sufficient thickness to contain the pressure loads of the helium gas 

and heat pipe fluid. The window must contain the helium pressure of 15 torr. 

If the distance, L « 20 cm (eq. (Cl)), then the maximum shear stress (p L/t ) 

gw 

will be 4 n/mm 2 (600 psl). This Is well within the yield stress for NaCl with 

Impurities added to Increase the yield strength (Table I). The potassium heat 

2 

pipe of reference 14 operates at a pressure of 4124 n/m (31 torr). Therefore, 

the pressure load on the titanium Is p = 31 - 15 = 16 torr, which for heat 

H 

pipe length (ref. 14), L' = 5.15 m, produces a maximum shear stress (p^L'/tp) 
of l.lxlO 2 n/mm 2 (1.6xl0 4 psl) well within the yield stress ( 1 . 55x1 0 5 psl). 

Using the materials, wall thicknesses, and dimensions already described, 
the specific mass ratio, a j/ a Q . can be calculated using equation (C14). 



= 0.25 
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Thus the reduction In wall thickness made possible by the window meteoroid 
bumper results In a 75 percent reduction In the specific mass. It should be 
remembered that the specific mass Includes only the radiating area. It does 
not Include the heat pipe working fluid nor other supporting structure. These 
parts will be Included when M t /M q (eq. (10)) Is calculated. 


In order to calculate 


Vo 


the ratios 


a'/a. 


and 


'/a 


/L-r-r %*«■ . 

T T oo 

tion (10) must be known. From the data of reference 14, a'/a = 1.26. 

0 0 

that the auxiliary material specific masses are the same (aj = a ), then 


In equa- 

Assumlng 


a'/a 

T T 


a'/a (a /a T ) 

o o o T 


Using this result plus a /a = 0.25, A„./A = 1.03, and 

To To 


a'/a = 1.26 In equation (10) yields M_/M = 0.69 
0 0 To 


Thus the 6PR results In a 
31 percent mass reduction from the comparable heat pipe radiator. The masses 
and areas obtained from these calculations are shown In Table II. 


For the hypothetical GPR the wall thicknesses, t w and t p were not 

reduced to the minimum required for meteoroid protection (t y + t = t Q /6.9). 

If t , t , t ■+ 0 then a T -» o and the minimum possible mass ratio M T /M 
w p s T K To 

will be obtained. From equation (10) the following Is derived. 


1 


a-j- -» 0 


1 + 


_o 

. i 


assuming 


a = a. 
0 


'. Therefore, for a /a = 1.26 and 
T oo 


A t /A q = 1.03, 


= 0.57 


-* 0 


H t /H = 0.69 obtained for 
T o 


Thus 

possible mass ratio 


V 


h t /m 

T 0 a T 


t : 
P 

= 0.57. 


t = 0.1 mm Is close to the minimum 
Therefore, little will be gained by 


reducing the wall thicknesses further. 


V. CONCLUSION 

This study was directed at predicting the performance of a new space radi- 
ator concept, the gas particle radiator (GPR). The GPR uses a gas containing 
emitting, submicron particles as the radiating media. The gas particle mix- 
ture Is contained between the radiator emitting surface and a transparent win- 
dow. There are two major advantages the GPR has over conventional heat pipe 
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or pumped loop radiators. First of all, high emlsslvlty Is achieved without 
the use of emissive coatings. Secondly, the GPR potentially has a lower mass. 

A radiation heat transfer analysis of the gas particle mixture yielded an 

_4 

expression for the GPR emlsslvlty, Cj. For a modest volume fraction (<p > 10 ) 

of submicron particles and gas thickness (D * 1 cm) It was found that the 

emlsslvlty was determined by the window transmittance. Thus the window becomes 
a critical element In the GPR concept. The window must have high transmittance 
In the Infrared and be strong enough to contain the gas particle mixture. 
Several candidate window materials are presented In Table I. The listed mate- 
rials are alkali halides and oxides of silicon, aluminum and magnesium. 

Another critical Issue Is maintaining a uniform particle distribution In the 
gas so that high emlsslvlty Is achieved. 

Besides determining the emlsslvlty, the window Is the critical element 
for making possible the lower mass for the GPR. The window acts as a "bumper" 
to provide meteoroid protection for the radiator wall. Thus the combined 
thickness (window + radiator wall) can be significantly less than the radiator 
wall thickness alone and still provide the same meteoroid protection. There- 
fore, the GPR can have a lower mass than a tube type radiator. 

The GPR was compared to a proposed titanium wall, potassium heat pipe 
radiator. For both radiators operating at a power level of 1.01 M w at 775 K It 
was found that the GPR mass was 31 percent lower than the heat pipe radiator. 

There are many design Issues that will have to be addressed and solved to 
make the GPR a viable space radiator. Two of the most Important Issues are; 

(1) Providing a seal between the window and the emitting radiator surface 

(2) A uniform and compatible thermal expansion of the window material and 
the emitting surface. 

Results of this study Indicate that the GPR can have a lower mass than 
tube type radiators without the use of emissive coatings. The window Is the 
critical element for the GPR and must provide high Infrared transmittance and 
sufficient structural strength. Also, to obtain the calculated high emlsslvlty 
a uniform particle distribution In the gas must be maintained. 


12 



APPENDIX A - EMISSIVITY OF GAS-PARTICLE RADIATOR 


Approximate the gas particle radiator as an Infinite flat plate with gas 
between the plate and a cover window. Figure (Al) Is a schematic of this con- 
figuration. In deriving an expression for the spectral emlsslvlty, e^, of the 
gas-particle radiator the following approximations are used. 

(1) Window at uniform temperature, T , and behaves In a diffuse manner 

w 

(2) Plate at uniform temperature, T^ and behaves In a diffuse manner 

(3) Gas at uniform temperature, T^ (no conduction or convection) 

(4) Negligible scattering of radiation by particles 

(5) Absorption coefficient, k^, of gas-particle medium Is uniform and Is 
determined by particle optical properties (gas Is transparent to 
radiation) 

(6) Steady state conditions, a/at = 0 

For the one-dimensional geometry shown In equation (Al) and neglecting conduc- 
tion and convection the energy equation Is the following, 


dx 


= 0 


(Al) 


where Qr Is the total radiative heat flux, (W/m^) 

Q r =/ o q x dX 


(A2) 


And q^ Is the specular heat flux (W/m m) . As a result of equation (Al) 


/•co /•oo 

, -/ < » •/„ < 


dx = constant 


(A3) 


n W 

where q^ Is the heat supplied to the plate and q^ Is the heat leaving the 
window. The hemispherical spectral emlsslvlty Is defined as. 


w 

c x ( V x) = — p 

p 

where Is the black body specific Intensity, 


(A4) 


B?(T ,X) = — 71T7T!t 
X p 5 1 hc/XkTp 


X V e 
13 


( hc/XkT \ 

(e V-y) 


(A5) 



Note that e Is defined In terms of the plate temperature, T . Also, the 
X P 

total hemispherical emlsslvlty Is the following. 


MV 



(A6) 


w 

In order to determine and Cj the heat flux leaving the window, q^ must 
be found. To obtain q? the radiative transfer equation must be solved with 
appropriate boundary conditions. First consider the boundary conditions. 
Taking a heat balance for a unit area of the window the following Is obtained 
(eq. (Al)). 


w 


w 

- + d. 


w 


w w 
+ 2c7nuV 

X \ 


(A7) 


w w 

Where q Is the heat from the gas Incident on the window, q Is the heat 

X 1 r 

reflected from the window back Into the gas, qlf the heat transmitted through 

x t 

the window and 2c^irB^ Is the heat emitted from the window to the gas and to 
the outside. In obtaining equation (A7) It has been assumed that no heat Is 
being lost by conduction and that no radiation Is Incident on the window from 
the outside. The total heat that leaves the window Is the following. 

<A8> 


Using equation (A7) to obtain 
yields the following, 


q^ and then substituting In equation (A8) 
X t 




1 


w 

<’ - o“> 


w B w 

v\ 


( A9) 


where p? = q. /q. Is the window reflectivity. Also, since It Is assumed 

X Xp x^ 

that the window behaves In a diffuse manner, the specific radiation Intensity 
leaving the window and entering the gas, I* , can be written In terms of the 

w x 

heat flux q" as follows. 

X 


w C T W „ . . T W W W ^ W D W /Ain\ 

q = / I cos 9 du : irl = P.d. + e.ifB ( Al 0 ) 

X f X XX X* X X 

0 / 0 0 0 1 
*'u=2ir 
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Now consider an energy balance for a unit area of the emitting plate, (eq. (Al)) 

(All) 


p p p 

\ - q x - \ 
A X Q X 1 


where q£ Q Is the heat leaving the plate and entering the gas and q^ Is the 
heat leaving the gas and striking the plate. Since the plate Is assumed to 
behave In a diffuse manner then a relation similar to equation (A10) for the 
window Is obtained. 


^ = f if cos © du = irlf 1 s p^q!^ + c f ir Bf 

0 0 x 0 X x 1 X X 


(A12) 


The boundary conditions (eqs. (A9) to ( Al 2) ) will be applied to the radiative 
transfer equation (ref. 1). Neglecting scattering, the radiation specific 
Intensity, I (S), at the window resulting from radiation within the solid 


1 

angle, do, (fig. (Al)) Is the following 

w o "M s > 

i” (S) = ij e X 

1 o 

Where Is the optical depth. 


r k (s) 

V. B ? exl> [-( K x-<) d <] 


-f 


(Al 3) 


V s )-/ k ? ds ’ 

0 


( Al 4) 


and Is the gas-particle absorption coefficient. Also, 


S = 


cos e 


(Al 5) 


Therefore, since T g Is assumed uniform (as a result B 9 Is uniform) and k 9 
is uniform, equation (A13) becomes the following. 

w d “ k f D / cos © a F -k 9 D/cos el 

r x - ij e x/ ♦ B 9 Ll - e x/ I < A16 > 

X 1 \ x 

The heat flux at the window resulting from radiation over all solid angles seen 
by a point on the window Is, 




cos e do 


( Al 7) 
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Substituting equation (A16) In (AT 7) yields the following 


w tP - „9- 

T x * 


(A18) 


where Is the gas transmittance and <* x Is the gas absorptivity. 


/ 

*u=2ir 




»oo q 

,2 I " k X S dS 


S=D 


e " cos 6 du = 2 D‘ / e * ^ = 2 E 3 (k q o) 


(A19) 


.*] 

o)=2ir 


cos 0 dco = 1 - t. =1 - 2E. 




(A20) 


and E (x) Is the exponential Integral. 

J 


E 3 (x) = 


00 

/■. 


e“ xu du 


( A21 ) 


w 

Similar to Equation (A18) for the window heat flux, q , the following result 

X 1 

Is obtained for the plate heat flux. 


P ,w - „g- 

q \, = * l \ T X * * B x°X 
1 o 


(A22) 


w 

If the boundary conditions, equations (A10) and (A12) are substituted for I In 

0 

equation (A22) and I q In equation (A18) then two equations for qj* and 

w 0 1 w 

q^ are obtained. These can be solved to obtain the following result for q^ . 

*1 *1 


\ = i W [ B ’(' - ; 0(’ * vJ * * 8 

1 1 ' P \ P \ T \ L 


W W P-2 1 

VWxJ 


( A23) 


Equation (A23) can now be used In equation (A9) to obtain the heat flux leaving 
w 

the window, q . Therefore, the spectral emlsslvlty (eq. ( A4) ) Is the following. 

K 


w 


1 - p. 


w 


c x = 


irB x 1 - 'x'Vx 


h JD - - *?>] 


B; 


c w b w 

X X 


r p i « w rP 

L B x 1 - p \ B x 


+ c x*x 


' B x 

1 - n 

B ? 


(A24) 
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In obtaining equation (A24) Klrchhoff's Law (ref. 1), 


g , and the relations 


a 


X 






(A25) 


w w w w w w 

p x + a x + T x = p x + c x + T x = 1 


(A26) 


were used, where p Is reflectivity, a Is absorptivity and t Is trans- 

X x x 

mlttance. 

The spectral emlsslvlty, c^, was derived under the assumptions given at 

the beginning of this appendix. Now make the additional assumptions that 

w w 

T = T = T and that the window has negligible emlsslvlty (t = 1 - p ). 
w g p x K 

Therefore, equation (A24) yields the following. 


w 


e. = 


1 - 


w p 

Wx 


72 [' - ; x(’ - ‘x)] 


(A27) 


To obtain the total emlsslvlty, c T , equation (A27) Is substituted In equation 
(A6) and the Integration Is performed. In order to simplify the g_ calcula- 

n y n_ 2 

tlon the following additional approximations are made * c x << 1 and 
« 1. Using these approximations will result In a value for lower then 

the actual value. Therefore, the total emlsslvlty, e^, will be a conservative 
estimate. 



Assumes ; 

•!«’. vM ^ 1 



(A28) 


Now consider the total emlsslvlty, g^.. As discussed In the Emlsslvlty of GPR 
section, there are two cases to be considered for the absorption coefficient, 
k^. For small particles (r d /X < 1) the absorption coefficient Is approximated 
as 



(A29) 


where K Is a constant and <p Is the volume fraction of particles In the gas. 
For. large particles (r./X > 1) It Is assumed that the absorption cross section, 

a O 0 

= irr^, so that, 
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( A30) 


9 _ 1 52 _ 
X " 4 r 


In order to calculate e T , an expression for the window transmission, x* must 

w 

be obtained. For these calculations the following result for x will be used. 


T. = 

T 

for X. < X < X 

X 

w 

SL - ~ U 

= 0 

for 

x < x ft and X > 


Also, x^ (eq. (A19), can be approximated as follows. (Ref. 1). 


( A31 ) 


\ = 2E 3( k X°) “«p(-'-< D ) (A32) 

If equations (A28), (A31), and (A32) are now substituted In equation (A6) the 
following Is obtained 


Cj - 


ITT 


W 





B X d) - 


( A33) 


For the case of large particles (eq. ( A30) ) Is constant and e T becomes 

the following. 


: T - ’w[’ - e "“][ F o- 


Vp 


- F. 


" X St T p 


6 


2.7 


(pD 


r 


d 



(A34) 
(A3 5) 


and F Is a tabulated function appearing In many radiation heat transfer 

0 — X I 

texts (ref. 1). 


F 


o-XT 



B x< x ’> 


d(X'T) 



dx 

x 5( e hc/kx _ ^ 


( A36) 


For window materials to be considered the largest value of X « 1 »im, and the 

1 -2 

highest temperature, T « 1200 K, It Is found that F Q .^qq vm k < * T ^ere- 

fore, the last term In equation (A34) can be neglected and 
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( A37) 


t ■ T « f o-x„t <’ - e ' B > • r > 1 

u p 

Now consider the case for small particles, (r./X < 1). Substituting equation 

d 

(A29) In ( A33) yields the following. 


‘ T ■ M [ F °-Vp ' Wp’ 


/ T X 
P U 


(- 4 * *>) 


T. x. <*V e *P Hf 


Pi 


ilr - 1 ' 


<l(*Tp)} 


(A38) 

Again the F Q T term can be neglected. In order to simplify the Integration 

the approximation exp (hc/XkT) » 1 Is made. This will be a good approxlma- 

4 

tlon as long as XT < 1.4x10 vm K (ref. 4). With this approximation the 
Integral can be written In terms of the Chi-Square probability function 
(ref. 15). Therefore, 


C T = T wV o-X. T 


90 

4 


up ir 


(r 1 l) Q [x B kT p (T * 1 


< 1 


(A39) 


where, Q(x |v) Is the Chi-Square probability function (ref. 15), 


Q(x‘ 


1/2 x 2 


dt 


( A40) 


r(n) Is the Gamma function and r(n) = (n - 1)! for n an Integer. The 
parameter y appearing In equation (A39) Is the following. 


3.6KkT 


T = 


he 


<pD 


(A41) 


In equation (A39) the term Q[2hc/X kT (y + 1 ) | 8] has been neglected since for 

a p _? 

than in 



APPENDIX B - TEMPERATURE DIFFERENCE BETWEEN PARTICLES AND SURROUNDING GAS 


Consider a particle to be approximated by a sphere of radius, r^. The 
surface of the sphere Is at temperature, T and Is surrounded by a gas at 


temperature, T , and moving with velocity, V. 
9 



rad 


For steady state conditions the heat transfer to the sphere will equal the heat 
radiated by the sphere. Therefore, 


h \' T g - T rad> " «‘sM T rad 


( T Jad - T l) 


(Bl) 


where h Is the heat transfer coefficient, A , Is the surface area of the 

8 2 4 

sphere, a(= 5.7x10“ watts/m K ) Is the Stefan-Boltzmann constant, Is the 
sphere emlsslvlty, and T Is the radiation sink temperature. Absorption by 

oo 

the surrounding gas, which will lessen the temperature difference, T^ - T^, 
is being neglected In equation (Bl). Since the largest temperature difference 
Is of Interest here, neglecting gas absorption Is a conservative approxima- 
tion. From Eq. (Bl) the temperature difference Is 


4T = < T g - T rad> ' 9 ( T Jad ' T l) 


(B2) 


The parameter, 


cfe 


6 = 


(B3) 


will determine AT. 

The heat transfer coefficient, h, depends upon the flow conditions. For 

the case of very small particles the mean free path of the gas molecules, 8. g , 

will be the same order of magnitude as r^. In this case the heat transfer to 

the particle will be described by free molecular flow. In the case of large 

particles, 8. « r^, and the heat transfer will be determined by continuum 

flow. First consider the free molecular flow case, (8. /r > 1). The energy 

9 d 

transfer to the sphere Is the following. 

Q = f (E 1 - E r ) dS ( B4) 
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Where Is the energy of the gas molecules striking the sphere and E r Is 
the energy flux of the gas molecules reflected from the sphere. From the 
definition of the thermal accommodation coefficient (ref. 16), a. 



(B5) 


where E Is the energy flux of gas molecules that have completely accommo- 
w 

dated to the sphere wall temperature, T Using equation (B5) In (B4) yields 
the following. 

Q = a J (E 1 - E w ) dS ( B6) 


Now consider expressions for E. and E . The energy fluxes, E. and E 

1 w 1 w 

consist of two parts; the kinetic energy flux, E', and the Internal energy, E", 
carried by the gas molecules. The kinetic energy flux, E 1 , Is given by 


/ CD r CO CO 

dC *i dC v/ 2 


V C X F dC 2 


(B7) 


2 2 2 2 

where c * c + c + c Is the gas molecule velocity squared and F Is the 
x y z 

distribution function. Assume a Maxwellian distribution function and since 


the flow velocity, V, Is small compared to the thermal speed, c 


m’ 


2 2kT 
c m “ m 

9 

the following results for F. 



Using equation (B9) In (B7) yields the following. 

E' = h(2kT) 

Where rt Is the gas molecular flux, 



( B8) 


(B9) 


(BIO) 


(Bll) 
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and n Is the gas number density. The Internal energy flux, E", Is the 
following. 


E" 



(B12) 


where J Is the number of degrees of freedom of the gas molecules which can be 
written In terms of the specific heat ratio, y* as follows. 

3 - Trf 1 (Bi3) 

Using equations (BIO), (Bll), and (B12) yields the following (ref. 16). 


E 


1 



+ E 


.. _ XJL-I— 
1 ■ 2(y - 1) 


"l kT g 


(B14) 


E 


w 



+ E 


« _ X_LJ 

w " 2(y - 1) 


"r kT rad 


(B15) 


Since the flux of Incident molecules, n^, must equal the flux of reflected 
molecules, n r , the following result Is obtained. 


E p - T - - 
t 1 “ "w ~ 2(Y - 1) "1 


6, 


*1 k < T g - T rad> 

(B16) 

/Birr 

V m 

(B17) 

▼ irm q 


Assuming T . Is uniform on the particle surface, the following Is obtained 
rad 

for Q by substituting equation (B16) In (B6). 


Q = 


g^x 


8( Y 


^1) "aV: 


8kT 

,m. 2 k < T , - T rad> A S 


(BIB) 


From this the heat transfer coefficient Is the following. 


h 


t 1) 

f = 8 (y - 1) 



(B19) 


Now consider the continuum flow heat transfer case. For low speed flow (Oseen 
Flow) the Nusselt number for a sphere Is the following (ref. 17), 



( B20) 
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for small Reynolds number. Appearing In equation (B20) are the continuum heat 
transfer coefficient, h c , and the gas thermal conductivity, \ t ^. From 
equation (B20) the following Is obtained. 



(B21 ) 


The thermal conductivity Is the following (ref. 18), 



(B22) 


2 

where Xg Is the velocity averaged cross section = 8ird g for hard spheres, 

where d Is the gas molecule diameter). If equation (B22) Is used In (B21) 

g 

and the result combined with equation (B19), the following Is obtained. 


lie _ 75 ( y - 1) J* ( ia\ 
h f " «( Y + 1) T8 \rj 


Where fi. 

9 


Is the gas mean free path. 


I 


1 



(B23) 


(B24) 


Now consider an estimate of the parameter, e, for the gas particle radiator. 

For helium at T = 1000 K, and a pressure, p = 0.1 atm, (n = p /kT ) and 
9 9 9 9 9 

Y = 5/3, equation (B19) gives the following result. 


hf = 1.2x10* a watts/m 2 K 

Also, using equation (23) the following Is obtained. 

h = 1 . 4x1 0 5 ( SL /r.) watts/m 2 K 
c g d 


Therefore, 


ac i o c 

6, = T- 1 = 4.8x10 — 1/K 

f hf a 


0 « r- 1 = 4.2X10" 13 

C h c 


( r d /l g ) 1 


1/K 


(B25) 


(B26) 


(B27a) 

(B27b) 
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Assume e $ /a * 1, also for p g = 0.1 atm and £ g * 10 -15 cm 2 , l g * 10 v m so 
that for r d * 100 ym, r d /a g * 10. Therefore, e f * 4.8xl0 -12 , and assuming 

e $ « 1 . e c 85 4.1xl0 -12 . Neglecting T* in equation (B2) yields the 
following. 


AT 


T 


rad 


6T 


3 

rad 


Using the above values for e f and e c the following Is obtained for 

T . = 1000 K L 

rad 


(B28) 


— I < 0.01 

rad 



< 0.01 


(B29) 


Therefore, for both free molecular flow and continuum flow the temperature 
difference will be less than 1 percent and can be neglected. 
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APPENDIX C - SPECIFIC MASS COMPARISON OF GAS PARTICLE AND TUBE RADIATORS 

Consider one segment of the gas particle radiator to consist of a metal 
tube surrounded by a transparent window as shown In equation (Cl). The window 
is attached to the tube by vertical members that are spaced a distance L 
apart. Therefore, the total mass of the radiator section, m^, not Including 
the mass of the radiator working fluid (wick and fluid mass for a heat pipe) 


Is the following 


1= p V + p V + v/p + <pp \ + p v 
T ww pp g^g dj s < 


Where p Is the density of the window material, V Is the volume of window 
w w 

material, p^ Is the density of the tube material, Is the volume of tube 

material, V Is the gas volume, p Is the gas density, <p Is the volume 
9 9 

fraction of particles, p^ Is the particle density, V $ Is the volume of the 
attachment member, and p^ Is the density of the attachment member. Referring 
to figure (Cl), the volumes are the following. 

V = irLt (2R + t ) ( C2) 

w WWW 

V = irLt (2R + t ) (C3) 

P P P P 


V* L 


K- 


2t R - t 
P P P 


V = irt 
s s 


(r 2 - R 2 ) 
V w p / 


Where t denotes material thickness, R denotes radius, and L Is the length 
of the radiator segment. 

Now define the specific mass, a^. In terms of the tube surface area. 


a T = 2^rL 

Substitute equations (Cl) to (C5) In (C6). Therefore, making the approxima- 
tion t « R and t « R the following result Is obtained, 
w w p p 

“T ■ 'vA© * Vp * * ~d) (’ 'I'l ( * <) t “ 

where 


D = R - R 
w p 
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For a simple tube radiator of thickness, t , and density, p , the specific 

0 0 

mass, a , will Include only the second term In equation (C7). 


m Q 

°o 5 2irR Q L = p o t o 


(C9) 


Now divide equation (C7) by (C9) so that the following result Is obtained. 


where, 


_ fa) + /D_\ fa) fa fa 

a o ** w V 0/ Pg Vo) lis \ t o/ p o VoJ 


(CIO) 


_ p w R w 

0 p 


_ p g * V 

W 9 " 2 Pft 


R 

11 + R* 

V 


(Cll) 


(Cl 2) 


_ 1 's/D\(, + V 

V s = 2^(l)(Wr- 


( Cl 3) 


for the case where the same tube material Is used for both the gas particle and 
simple tube radiator (p Q = p p = p $ ) then equation (CIO) becomes the following. 


fl (jd\ ( D_\ . fa Ji 

“o = "AO Vs W V 




(C14) 


Where, since p = p 
r s r o 


= 2® * R p ; 


(Cl 5) 


If a flat plate radiator Is being considered rather then a cylindrical 
geometry then, 


- w 
Ww = 'o 

( Cl 6) 

p o \ p o / 

( Cl 7) 
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(C18) 



D = distance between emitting plate and window 
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TABLE I. - CANDIDATE WINDOW MATERIALS 


Material 

Transmission range 

Transmittance 

Material 

thickness, 

mm 

Density, 

gm/cm 3 

Melting 

point, 

•c 

Water 

solubility, 
gm/100 gm 

(Type) 

Yield stress. 
N/mm 2 (lb/in 2 ) 

x £, 9 pm 

x u , um 

Fused 








silica 

0.25 

3.5 

>0.9 

2.65 

1610 

Insol . 

— 

(Si0 2 ) 



(5) 





Sodium 







(Compressive)^ 

chloride 

.35 

15 

>.9 

2.165 

801 

35.7 

0.69 to 41.3 

(NaCl) 



(5) 




(100 to 6000) 

Silver 








chloride 

.4 

20 

.8 

5.56 

455 

0.0021 


(AgCl ) 



(.5) 





Calcium 








fluoride 

.25 

8 

.95 

3.18 

1423 

.0017 


(CaF 2 ) 



(1-id 





Lithium 







(Shear)l 2 

f 1 uori de 

.2 

5.5 

.95 

2.635 

845 

.27 

1.2 to 6.0 

(LiF) 



0-3) 




(174 to 870) 

Sodium 

.3 

10 

>.9 

2.558 

993 

4.22 

- - 

fluoride 



(2.16) 





(NaF) 








Potassium 







(Compressive)^ 

bromi de 

.5 

20 

>.9 

2.75 

734 

102 

0.69 to 40 

(KBr) 



(4) 




(10 to 5800) 

Cesium 







(Compressive)^-^ 

bromide 

1 

25 

.9 

3.04 

636 

123 

22.5 

(CsBr ) 



(10) 




(3270) 

Aluminum 








oxide 

.3 

4 

.85 

4.0 

2316 

Insol . 


(AI2O3. 



(0.5) 





crystal ) 








Magnesium 







(Shear)l 2 

oxide 

.6 

6 

.85 

3.77 

3223 

Insol . 

36 

(MgO, 



(1-9) 




(5225) 

crystal ) 








Cesium 








iodide 

1 

40 

.9 

4.51 

626 

160 


(Csl, 



(5-10) 





crystal ) 








Heavy-metal 








fluoride 








glasses 








(HMFG) 








HF 4 -BaF 2 - 

LaF 3 *AlF 3 

.3 

7 

.9 

5.88 

312 



BaFe -ZnFp* 








LaF 3 *TnF 4 

.3 

9 

.9 

6.2 

357 


— 













TABLE II. - COMPARISON BETWEEN HEAT PIPE RADIATOR AND GAS 
PARTICLE RADIATOR 


[Radiator conditions; radiated power = 1.01 MW, radiator 
temperature = 775 K]. 



Radiator 

area, 

A, n? 

Radiating 
surface mass, 
*Vad = “A, kg 

Auxiliary mass, 
M a = b'A, kg 

Total mass, 
M rad + M a> kg 

Heat pipe 
radiator 
(ref. 14) 

67 

170 

214 

384 

Gas 

particle 
radi ator 

69 

44 

221 

265 


Titanium wall potassium heat pipe data (ref. 14) : Emissivity, 


c 0 = 0.9; radiating wall thickness, t Q = 0.6 mm; area 
redundancy, r = 321/360 

Gas particle radiator data: Emissivity, cj = 0.84, NaCl window 
U u = 15 pm, t w = 0.9); wall thicknesses, t p = = t $ = 0.1 mm; 

area redundancy, r = 321/360; redundant area emissivity, 

= 0.3; helium pressure, p g = 15 torr. 

Calculated results used to obtain GPR results: 


a T 

— = 0.25 , 
a o 




= 1.03 



1.26 
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Figure 1. - Gas particle radiator concept. 
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Figure 2. - Total emissivity for small particles 
(r d < 1 pm, T p < 1200 K). 




Figure 3. - Total emissivity for small particles (r d <1 |Jm, 
Tp < 1200 K). 
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(a) Cross section of heat pipe radiator of 
Ref. 14. 
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(b) Heat pipe radiator modified to produce a gas-particle radiator. 

Figure 4. - Cross-section of heat pipe radiator of Ref. 14 and 
modification to produce a gas-particle radiator. 
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